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I. INTRODUCTION 



Recently, the central exclusive production of y c charmonia has attracted a lot of attention 
from both experimental [l|-[3| and theoretical |H, [H| sides. At the moment, such a process 
provides the unique opportunity to test the QCD diffractive Kaidalov-Khoze-Martin-Ryskin 
(KKMR) mechanism |6[ based on the /c r factorisation incorporating nonperturbative small- x 
gluon dynamics described by the unintegrated gluon distribution functions (UGDFs) against 
accessible data. 

The final hadronic central system is, however, rather complicated and composed of three 
spin states Xc(0 + , 1 + , 2 + ), and there is no yet any reliable way to measure them separately. 
Such a measurement would significantly reduce an overall theoretical uncertainty, as dif- 
ferent Xc spin contributions come from different phase space regions 0^9|- Indeed, in the 
asymptotical forward limit xA^ + , 2 + ) are strongly suppressed with respect to Xc(0 + ) due to 
the J z = selections rule However, corrections to this asymptotics turn out to be im- 
portant in the total cross section leading to a noticeable contribution of the Xc(l + , 2 + ) states 
to the observable signal from radiative Xc decays 0,0,11]. Generally, separate measurements 
of Xc states in diffractive production would impose more strict bounds on the production 
mechanism under consideration and extract some new information on the underlying QCD 
dynamics (El, [sj . 

One of the ways is to measure the characteristic differential distributions like meson p± 
distributions or angular correlations of the outgoing protons. However, such distributions are 
rather sensitive to UGDFs and cuts on kinematical variables 0,0], and existing experimental 
techniques do not allow to reconstruct such observables with sufficient precision. Another 
possible way is to look at observables related with meson polarisations and their decay 
distributions. 

The goal of the present paper is to analyze polarisation effects in the central exclusive 
production of Xc charmonia. Such effects can be potentially identified by measuring the 
angular distribution of J/ip mesons from radiative decays of Xc{J + ) giving more detailed 
information on the partial meson helicity contributions. Moreover, certain combinations of 
polarisation observables can be less sensitive to unknown nonperturbative effects leading to 
unique opportunities for model-independent analysis of diffractive processes. 

The paper is organized as follows. In Section II we give general expressions for the 
amplitude of the central exclusive Xc(l + , 2 + ) production with all necessary notations. Section 
III contains discussion of the rapidity dependence of the hard subprocess amplitudes for 
different Xc polarisations. In Sections IV and V we derive the angular dependence of outgoing 
J/ijj mesons in the helicity frame in terms of the diffractive Xc production density matrix 
p\y. Sections VI and VII are devoted to discussion of gluon off-shell effects and absorptive 
corrections. Section VIII contains the presentation of the main results, including differential 
distributions of polarised x c (l + ,2 + ) and angular correlations of the outgoing J/ip mesons. 
Finally, in Section IX we give a set of concluding remarks and present a discussion of the 
final results and theoretical uncertainties. 
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II. AMPLITUDES OF EXCLUSIVE P-WAVE CHARMONIA PRODUCTION 



According to the Kaidalov-Khoze-Martin-Ryskin approach (KKMR) p, we write the 
amplitude of the exclusive double diffractive color singlet production pp — > ppXcj as 

KA3*9* — S 2 ^cic 2 cv f ,2 t rc x c 2 fg,l( Xl ' X l ' ?0,i) Ql,t > *0 ( X 2, ^2 5 <?0,t > ?2,t) ^0 , . 

- 2 •* Af^T s i d ** v » 4j^[ t • C") 

where /°f ( 

^i? 9o,t> t> are the off-diagonal unintegrated gluon distributions for "active" 
gluons with momenta = XiPi + q ijt and color indices q, and the screening soft gluon with 
small fraction x' x,, g — 9o,t) ^ an d ^ are the s pi n an d helicity of a produced meson with 
momentum P = g x + q 2 and mass M in the center-of-mass frame of colliding protons and z 
axis directed along the meson momentum P, respectively. Vj\ C2 here is the hard g*g* — > Xcj 
production amplitude for charmonium with quantum numbers J, A. 

The structure of Vj in Eq. (12. ip (we omit here the color and polarisation indices for 
simplicity) is determined by gauge invariant amplitude Vj V for the off-shell gluon fusion 
process, and given in terms of its projection onto the gluon polarisation vectors as 0- 



^ = n>-VT = - — — K^, <£Vr = <£VJ*» = ■ ( 2 -2) 

In order to define the polarisation states of x c (l + , 2 + ) mesons, one has to fix a specific 
frame, in which those are uniquely identified and can be measured experimentally. The 
commonly used polarisation frames are the helicity, target, Gottfried- Jackson and Collins- 
Soper frames (see, e.g. [loll). In this work we primarily focus on polarisation effects in the 
helicity frame. 

According to Refs. 0, IsJ, in the considered frame we introduce the time- like basis vectors 
^1,2,3 satisfying n^n^g^ u = g a p (with Uq = P^/M) with collinear n 3 and P vectors (so, we 
have P = (E,0,0,P z ), P z = |P| > 0), and 

J _ n i n\ „P 1 



< = (0,1,0,0), n£ = (0, 0, 1, 0), 4 = — (|P|,0,0,.E), |P| = VE*-M\ (2.3) 

Then the hard production amplitudes for the axial-vector J = 1 (A = 0, ±1) and tensor 
J = 2 (A = 0, ±1, ±2) charmonia in the helicity frame read (for more details, see Refs. 0,0]) 



A ^ y V MvriV, |P t |(M 2 - g 2 , - g| t ) 2 I ^ 

■%i,t + 92,*) |[qi,t x <te,t] x m| sign(Q 2/ ) sign(cos^) 

(1 - l^l)(9i 2 t + ?L) |[qi,< x 02,*] x n 3 | sign(Q y )sign(sin^) 



i\M(<llt - 92,t)(?i,t?2,t)sign(siii^) + 
+ (2-4) 
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VjLlx - 2ig 2 5™ ] J M | Pt | 2(M2 [ ^ _ q 2^2 x ( 2 - 5 ) 

6M 2 i\\\(qj t - ql t ) sign(Q y )<! |[q M x q 2 , t ] x n x | (1 - |A|) sign(sin ifi) sign(cos ip) + 



2 |[qi, t x q^] x n 3 | (2 - |A|)} - \lq\ t q\ t + (ql t + g£, t ) 
|3M 2 (cos 2 ^ + 1)A(1 - |A|) + 6M£sin(2^)A(2- |A|) sign(sin^) sign(cos^) + 

v / 6(A^ 2 + 2E 2 )sin 2 ^(l- |A|)(2- |A|)' 



where ±Q y are the y-components of the gluon transverse momenta 51/2,* in considered co- 
ordinates, if) = [0 ... 7r] is the polar angle between P and the c.m.s. beam axis, and 



[qi,t x qa,t] x rixl = \jq\ t q\ t - {qx^t) 2 \cosip\, 
[qi,t x q 2 ,t] x n 3 | = jrJqJ^t - {qi,tQ2,t) 2 I sin ^l, 



M 

l p t| 2 = |qi,<| 2 + I<l2,t| 2 + 2|q ljt ||q 2it | cos 



where is the angle between fusing gluons. Amplitudes (I2.4p and (I2.5P explicitly obey gauge 
invariance and Bose symmetry properties with respect to the gluon momenta interchange 
qi -H- q 2 . 

Another important feature is that the amplitude of axial-vector charmonia production 
turns to zero for on-shell gluons, i.e. when q\ t = 0, q 2 t = due to the Landau- Yang theorem 
(see, e.g. Ref. Q). The gluon virtualities (transverse momenta) provide a leading effect in 
the diffractive x c (l + ) production, and, therefore, cannot be neglected j3, 0]. This is the 
striking difference between the fc r factorisation approach under consideration and collinear 
factorisation which forbids production of axial- vector states [3]. 

As was demonstrated in Refs. 0, [1] amplitudes for diffractive Xc(l + ,2 + ) turn to zero 
in the forward limit, i.e. when qi^ = — q^.t = Qo,t) due to symmetry relations. This is a 



direct consequence of J z = selection rule [11] saying that CEP of higher spins J = 1, 2 



is strongly suppressed in the forward limit. Like the gluon virtualities in the Xc(l + ) case, 
the off-forward corrections provide a leading effect in the case of Xc(l + ; 2 + ) production, and 
cannot be neglected in the integrated cross section. In particular, they lead to a substantial 
contribution of x c (2 + ) meson CEP Q|, close to that from Xc(0 + , 1 + ) mesons. 

Let us now turn to the discussion of the polarisation effects in diffractive 1 + , 2 + charmonia 
production and first start from analytic investigation of helicity amplitudes. 

III. RAPIDITY DEPENDENCE OF g*g* -> X c(J) AMPLITUDES SQUARED 

For the purpose of illustration 1 , it is interesting to look at the y-dependence of the 
9*9* Xc(J = 1,2) (hard) subprocess amplitudes (12.41) and (12. 5p for different meson helic- 



1 Please, note that the content of this Section aims to analytical illustration of specific rapidity dependence 
of the hard subprocess g*g* — ¥ Xc amplitudes squared, not the diffractive pp —> ppXc amplitude squared. 



4 



ities A = 0, ±1 and A = 0,±1,±2, respectively. It is convenient to express them in terms 
of the transverse 3-momenta of fusing off-shell gluons |qi,t| and |q2,t|, and the angle be- 
tween them in the center-of-mass frame of colliding nucleons with the z-axis fixed along 
meson momentum P. In this case, summing the matrix element squared |V^ J | 2 over meson 
polarisations A up to some constant normalization factor N J we get: 



J=i — *tJ=i 



S J=1 = N' 



S J=2 = N J=2 





2 


Q2,t 


2 




2 + 


qi,t 


2 ) 2 sin 2 + M 2 (|q lit 


2 + 




2 -2 




|qa,t 


COS</>) 


(|qi,* 


2 + 


02,* 


2 + M 2 ) 4 



|qi,* 


2 


Q2,t 


2 


3M 2 M 2 + (|q M 


2 COS + 


qi,* 


2 cos + 2 q lj( 


|q2,*l) 2 




2 + 


qj2,< 


2 + M 2 ) 4 



(3.1) 

These sums are proportional to the ones derived by Kniehl et al. in Ref. [l2j (see right below 
Eq. (27)). The distinction is onl y du e to different normalizations of gluon polarisation vectors 
e M = Q1/2 t/|qi/2,t| use d m Ref. uM an d light-cone vectors used in our calculation. The 
sums (13. ip are y- independent. Then, for partial polarisations of the x c (l + ) meson we have 

j rl2 = s j =1 AMlco^y 



v ( 71 /T 2 ^ 1--* ^ 



(M 2 cosh 2 y - M 2 ) (EM 2 + A) 

, v j=i ,2 ^ J=1 ^ 2 [®( M i y-M 2 )- a] 

1 A=±l1 2 (M 2 cosh 2 y-M 2 )(SM 2 + ^)' 1 ' ' 

where 

-4 = (|qi,t| 2 + |qi,t| 2 ) 2 sin 2 0, B = |q M | 2 + \q 2 ,t\ 2 - 2|qi,t||q2,t| cos<f). 
Analogously, for the tensor x c (2 + ) meson we have 

J=2 2 _ J=2 K[Mi cosh 4y + 4(M 2 + M 2 )M 2 cosh 2y + 3M 4 + 4M 2 M 2 + 2M 4 ] 
'^ A=0 ' = 5 Q(M 4 cosh4y + AM cosh2y + 3M 4 - 8M 2 |PJ 2 ) ' 

J= 2 _ S J=2 3MiM 2 [\~P ± \ 2 Ml cosh 4y + 4V(cosh2y+l)-\~P ± \ 2 Ml] 
I A=±il -— ^ Q(M| cosh4y + 4A4 cosh2?/ + 3M 4 - 8M 2 |PJ 2 ) ' '" K ' >) 

j=2 1 2 s J=2 3M * [ M l ( cosh % + !) - 4M 2 M 2 cosh 2y - 2P + 2M 4 ] 



A=±2l 



2 Q(M 4 cosh4y + 4A4cosh2y + 3M 4 -8M 2 |PJ 2 ) ' 
where for compactness we have introduced the following short-hand notations: 

TZ = (|qi it | 2 cos0+ |qi i4 | 2 cos0 + 2|q M ||q 2 ,t|) 2 , 

Q = 3MlM 2 + H, P=|Pj_| 4 -^, M = \P ± \ i -M\ 



For illustration in Fig. [TJ we show the helicity amplitudes squared as a function of me- 
son rapidity at some fixed kinematical point. Apparently, they obey a quite non-trivial 
strong ^/-dependence in the vicinity of y — which will result in the "peaked" structure 
(maxima/minima) of the differential cross sections da^/dy of double diffractive polarised 
Xc(l + ,2 + ) production (see Results section). 
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FIG. 1: Illustration of the meson rapidity dependence of the normalized subprocess g*g* — > Xc(l + ) 
(left panel) and g*g* — > Xc(2 + ) (right panel) amplitudes squared for different meson helicities A at 
arbitrarily fixed kinematical variables Iqi t| = |qi t| = 0.5 GeV and <fi = 0.2. 



It follows from Eqs. ( 13. 2 p and ( 13 .3p that in the forward limit, corresponding to qi t = — q 2j t, 
helicity amplitudes V^Tq and V^J7q ±1 turn to zero. So, the total signal is dominated only 
by maximal helicity contributions, i.e. by A = ±1 for % c (l + ) and A = ±2 for % c (2 + ). 
Therefore, nontrivial "peaked" rapidity dependence of the cross section around y = 0, as 
long as absolute values of non-maximal A contributions, can be served as a measure of 
non-forward corrections in the production process. 

Let us consider how to measure such polarisation effects in radiative decay channel 
X c (l + ,2 + ) — > J/ip + ry and understand how these effects manifest themselves in differen- 
tial distributions of outgoing J/ip meson. 



IV. DIFFRACTIVE X cj PRODUCTION DENSITY MATRIX 

The cross section for the unpolarised Xcj production in the 3-body reaction pp — > ppXcj, 
where J = 1,2, can be written as 

< = E ^A, <> = < " P J XX = Y 3 I d3PS • ^A'({-..}) > 

A=-J J 

where p xx , is the hermitian helicity density matrix p xx , = p J x , x . In the differential form it 
can be written in terms of invariant diffractive amplitude as 

d PM{~}) = ^T L =-±r-\ E ^a 1 a 2 ^ 1 a' 2 a({---})^a 1 a 2 ^ 1 a^({---}) ■ (4-2) 

Above {...} is an abbreviation for a four- dimensional phase space point. In our case, because 
of nucleon helicity conservation (Ai = \[, A 2 = A' 2 ), the nucleon helicities can be suppressed. 
In general, p xx , is a function of kinematical variables (xF,ti,t2, <fi pp ) (or any other equivalent 
set of variables). 

Invariance of the amplitude A4 under reflection in the production plane by the action of 
the operator exp(— mJy)V , where V is the parity operator, leads to the relation 

P^A,-A' = (-1) A - A 'P3['A- (4-3) 
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Thus, for example, in the case of J = 1 there are only four independent components in p xx r. 
Pq , p\ x , p\ , and p_n. The diagonal elements of the density matrix are real. 



V. ANGULAR DISTRIBUTION OF J/ip MESON FROM X c(J) DECAY 

Let us consider the central exclusive production process pp — > ppXcj followed by the radia- 
tive decay Xcj J/ip + 1- We shall limit ourselves here to the narrow- width approximation. 
Below we follow notations in Ref. [l3[ • Let 9 and <fi be the polar and azimuthal angles of the 
J/ijj meson in the Xc{J + ) rest frame (this is so-called helicity frame). Then, the differential 
cross section of the J/ip production in the sequential process pp — > ppixcj J/4> + 7) can 
be written as 

da J J 

-£± = Bj( XB j^J/ih)-W{0,<l>) } W(9,<P)= fixMxiW), (5.1) 

A,A'=-J 

where dfl = d cos 9d<p, Bj denotes the branching fraction of Xcj — > J/ipl, W(9,<f)) is the 
angular distribution of the J/ip meson, p AA , is the integrated density matrix corresponding to 
Xcj production process (14. ip . and A xx , refers to the Xcj decay process and allows to describe 
the decay angular distribution 

AM$ - * ] = — — ■ (5 - 2) 

where Ai = 0, ±1 and A 2 = ±1 are the helicities of the J/ip and 7 in the Xcj rest 
frame, respectively. The transition matrix element of the radiative decay process Xcj(A) — > 
J/ip(Xi,9,(j)) + / y(X 2 ,TT — 9,(f) + it) in the considered frame of reference is [l4j 



'2 7-1-1 

(X 1 ,X 2 ,e,<f>\T\J^) = \l— — T/ 1>A XW-^M), (5.3) 



lX^ = \/2^i^i.^,0,0|T|J,A>| 



47T 

A=Ai— A2 



Let r], 7/1, and r\ 2 (J, Ji, and J 2 ) are the parities (total angular momenta) of the Xcj, J/ip, and 
7, respectively. Angular- momentum conservation imposes the selection rule |Ai — A 2 | < J- 
Due to theparity conservation the decay amplitude T XiX2 satisfies the following symmetry 
property [l4| 

T_ J Alj _ A2 = VVM-l) Jl+J2 - J T{ lX2 = (-l) J T/ iA2 . (5.4) 

Thus, an independent set of T x Aa components for J = 1, 2 reads 

J = 1 : t\ = T x l = —T_ x _ x , t\ = Tq _ 1 = — Tq 1 ^ , (5.5) 

J = 2: t 2 = ^1=^1,-1' *i = T o,-i = T o,i' ^2=^-1 = ^1. (5.6) 
The matrix 

DJ m>m( a >P>l) = {j^ m '\ D (a,P,l)\j,m} = exp(-i-fm')d J mlm (/3) exp(-iam) (5.7) 
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is the representation of the rotation operator 

D(a, f3, 7) = exp(-i7J 2 ) exp(-i(3J y ) exp(-iaJ x ), 



(5.8) 



with a, (3 and 7 being the Euler angles, in the eigenstates \j,m) of J 2 and J z . The well- 
known c?-functions d 3 m , m (f3) = (j,m'\ exp(—i(3J y )\j,m) may be evaluated from the Wigner's 
formula [151 ] 



min(j+m ,j— ml) 



d j / (p) = J2 ( l) k - m + m ' V U + m V-U - m V-U + - m')\ 



fc=max(0,m— ml) 
2j+m— ml — 2k 



x I cos - 



Orthogonality condition 



sin — 
2 



k\(k — m + rn')\(j + m — k)\(j — m' — k)\ 

2k—m+m' 



Atx 



2j + 1 



provides the normalization of A J XX , 



(5.9) 



(5.10) 



AA'VT 



#AA' 



(5.11) 



so that, upon integration over the solid angle, Eq. (15. ip reduces to the narrow-width ap- 
proximation formula 



(5.12) 



Combining all above ingredients together, we get the angular distribution of the J / ip meson 
( 15. ip in the general form 



W 



j=i. 
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rl cos 2 9 + — sin 2 9 



rl sin 2 + -A(l + cos 2 9) 



V5sin(20) (rl-^) [Re(p} j0 ) cos0 - Im(pi ) sin, 



— sin 9 r n — 



.1 ' 1 



[Mp\,-i) cos(20) - bn(pLi) sin(20)] I, 



(5.13) 
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J =2 1 



W 



+ Pli 



+ P2.2 



5a 



J=2 



4:71 



\rl(3 cos 2 9 - l) 2 + Irl sin 2 9 cos 2 9 + ^r 2 sin 4 9 
4 u 2 8 



3r 2 sin 2 cos 2 9 + ^ (4 cos 4 0-3 cos 2 + 1) + ^ sin 2 0(cos 2 + 1) 

9 r 2 2 

-r 2 sin 4 9 + - 1 - sin 2 fl(cos 2 5 + 1) + (cos 4 + 6 cos 2 + l) 
4 2 8 v ' 



+ — sin 2#(2r 2 (l - 3 cos 2 9) + 2r\ cos 29 + sin 2 0) [Re(p 2 ) cos — Im(p 2 ) sin 0] 
- - sin 2 9(6rl cos 2 9 + r 2 (sin 2 9-3 cos 2 9) - r 2 sin 2 0) [Re(p 2 _ x ) cos 20 - Im(p 2 _ x ) sin 2 
H — sin 3 6>(6r 2 — 4r 2 + r 2 ) cos 9 [Re(p 2 _ x ) cos 30 — Im(p2 -1) srn 3 



+ - sin [Re(p2 _ 2 ) cos 40 — Im(p2_ 2 ) s i n 4 



+ sin 2 #(2r 2 (3 cos 2 - 1) - 4r 2 cos 2 9 + r 2 (cos 2 + 1)) [Re(p^ ) cos 20 - Im(p 2 ) sin 20] 



1 



sin 20 (6r 2 sin 2 9 + 4r 2 cos 2 9 — r 2 (cos 2 9 + 3)) [Re(p 2 1) cos — Im(p2 x ) sin 0] 



(5.14) 



where 



UJ\2 



V If J I 2 



(5.15) 



are positive numbers satisfying ^2 x =o r x = 1 to be determined experimentally. In Ref. [16 



these values were extracted from the Fermilab E835 data on fractional amplitudes of the 
electric dipole (El), magnetic quadrupole (M2), and electric octupole (E3) transitions in 
the exclusive reactions pp — > Xcj J/'i'l ~ e + e~7, with J = 1, 2. They are listed here 



0.498 ±0.032, 
0.075 ±0.029, 



0.502 ± 0.032, 
0.250 ± 0.048, 



0.674 ±0.052. 



(5.16) 
(5.17) 



The corresponding values for a pure El transitions read r\ = r\ = 0.5, Tq = 0.1, r\ = 0.3, 
and r\ = 0.6. Having these values we only need to calculate the independent components of 
p J xx , to determine the J/ip angular distribution W J (9,(p) completely. 

Analogously to Eq. (14. ip . we can represent the angular distribution function as 



W" 



1 

2s 



d 3 PS -dW J (9, 4>; {...}) 



(5.18) 



Such a representation will allow us to calculate correlations between kinematical variables of 
the central system as a whole (for example, rapidity y, momentum transfers ti i2 and angle 
between outgoing protons PP ) and angular variables of produced J/xjj meson. Distributions 
dW(9, 0; {•••}) contain all the information about diffract ive Xc production mechanism, meson 
polarisations and decay process. 
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VI. GLUON OFF-SHELL EFFECTS 



Let us comment now on the importance of gluon virtualities (transverse momenta) in 
the hard subprocess vertex g*g* — > Xc{J + ) implied by the genuine /c t -factorisation scheme. 
This question was originally investigated in Ref. (9[ when considering the central exclusive 
production of Xc(0 + ) meson. In particular, it was found that the on-shell approximation 
(M^ c 3> q\ 2 ) relying on the neglecting terms of the order of q\ 2 /M^ is too crude for the 
production of low mass systems with mass being of the order of a few GeV. In the Xc(0 + ) 
CEP case, inclusion of the gluon virtualities q\ 2 m the hard matrix element leads to a 
noticeable decrease of the total (bare) CEP cross section by a factor of 2 - 4 depending on 
the model of UGDF used. 

This off-shell effect is even more pronounced in the case of the axial- vector % c (l + ) pro- 
duction, as it is dominated mostly by relatively large gluon virtualities and vanishes in the 
strict on-shell case. Even when included in the numerator of the vertex (see, e.g., Eq. (12.41) ) 
the "small" terms q\ 2 are sometimes neglected in its denominator. Such an approximation 
M£ c ^> q\ 2 increases the total Xc(l + ) CEP cross section by almost a factor of 7 in absolute 
normalization making a large effect in the Xc(l + )/Xc(2 + ) ratio in the observable radiative 
decay channel [§[ (see also Table [I] below) . This resulted in a large discrepancy (in both nor- 
malisation of the total cross sections and relative spin contributions into observable signal) 
between our previous results in Refs. @-|9j], where the exact hard matrix elements with off- 
shell gluons were used, and HKRS results in Refs. 0, [H| , where gluon on-shell approximation 
M 2 ^ 3> q\ 2 was adopted 2 . Unlike the absolute normalisation of the cross section, the shapes 
of the differential distributions are not strongly affected by such an approximation. 

One should note here that the description of the Tevatron data on the central exclusive 
Xc production with updated results, including off-shell effects, becomes more difficult as the 
total cross section in the radiative decay channel Xc(J + ) —> + 7 gets underestimated by 
about a factor of two with respect to the Tevatron data [8], if one uses values for the gap 
survival factors calculated in Refs. 0, @] • However, such results are still within an overall 
theoretical uncertainty from the soft physics involved, UGDF models in the QCD mechanism 
under consideration [8j, and also from unknown NLO corrections to the g*g* — > Xc(l + ) 
vertex, which are implicitly, but only partly, taken into account by the /^-factorisation 
approach used. 



VII. SOFT RESCATTERING CORRECTIONS 

In Refs. for simplicity the whole analysis of observables has been performed adopting 
an approximation when all soft rescattering effects are effectively absorbed into the effective 
gap survival factor (S^ s ) taken into account multiplicatively. However, it is known from 
earlier analysis of the central exclusive Higgs production jsj (and later for unpolarized Xc 
production Q) that the soft rescattering may change angular distributions crucially, so it 
is not enough to include the gap survival factor as a constant normalization of the cross 
section. 

In order to get some insight on the soft rescattering effects and include them dynamically, 



2 We are grateful to L. Harland-Lang for the significant help in understanding the major part of discrepancies 
between our calculations. 
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in the first approximation we use the simple one-channel model for the eikonal part of the 
survival probability calculated as 



cik(Pllt ' P2,t) " \M^{p^p 2 ,r ( } 



where Pi/2,t are the transverse momenta of the final protons, M bare {pi,t-,P2,t) is the "bare" 
amplitude for the Xc CEP, defined in Eq. (12.11) . and M res (pi,t,P2,t) is the rescattering 
correction which can be written in the form (see, e.g., Eq. (13) in Ref. [l7| ) 



where ImM (s) = sa^{s) (the real part is small in the high energy limit), B is the t-slope 
of the elastic pp differential cross section, and b ~ 4 GeV -2 is the t-slope of the proton form 
factor. At the Tevatron energy we have 0"*°* ~ 80 mb and B ~ 17GeV~ 2 . 

In Table [I] we have collected results for differential cross section da Xc /dy(y = 0) for 
different spin states, as well as their contributions to the radiative J ftp + 7 decay channel 
d<7j/^ y /dy(y = 0) at Tevatron. For illustration, here we use only the perturbatively modeled 
KMR UGDF 0,0, which include the Sudakov form factor and GRV94HO gluon PDF (20 



The results are presented for both versions of the absorptive corrections - by using the one- 
channel approximation according to the eikonal model 3 (17. ip . (17. 2p . and in the factorized 
form with p r averaged (effective) gap survival factors (S 2 S ) obtained recently in Ref. 

(S c 2 ff ( Xc (0 + ))> ~ 0.0284, (S 2 ff ( Xc (l+))) ^ 0.0735, (5 2 ff ( Xc (2 + ))) ~ 0.0539 , (7.3) 

where it is assumed that the enhanced suppression factor (S^ nh ) = 0.49 is the same for all 
spin states. 

With the on-shell approximation M 2 c 3> qf 2 we basically reproduce the HKRS observable 
signal ~ 0.7 nb [5]. However, we see from Table U that keeping the exact matrix elements 
for the hard subprocess g*g* — > Xc(J + ) (without the approximation M 2 3> q\ 2 adopted in 
Refs. @,[5[) it is rather difficult to get the Tevatron experimental value j2( da exp / dy\ y= o(pp — > 
pp(J/ip + 7)) ~ (0.97 ± 0.26) nb, even by moving lower cut on the transverse momentum of 
the gluons in the loop (q^f ) 2 down to the minimal perturbative scale 0.36 GeV 2 in GRV94HO 
gluon PDF [20]. This practically means that the central exclusive charmonia production is 
very sensitive to and seemingly dominated by yet unknown non-perturbative gluon dynamics 
at small x and q t . 

However, if we exclude the enhanced-diagram correction factor S cn h, which stands for 
QCD factorisation breaking in the considered QCD diffractive mechanism, we basically 
reproduce the Tevatron data and get the observable signal ~ 0.94 nb compared to the 
experimental value 0.97 nb. This may provide a clue that the QCD factorisation breaking 
effects may be overestimated and need further more careful attention. 

One-channel eikonal model (I7.2p is expected to provide the main effect on the shape of 
the distributions in the azimuthal angle $ between outgoing protons shown in Fig. [2j Other 
distributions in meson rapidity y and momentum transfers along each proton line £x,2 are 



3 As for enhanced-diagram correction factor S^h, we are grateful to L. Harland-Lang for providing us with 
a grid for its numerical calculation as a function of transverse momenta. 
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TABLE I: Differential cross section da Xc /dy(y = 0) (in nb) of the exclusive diffractive produc- 
tion of Xc(0 + ,l + ,2 + ) mesons and their partial and total signal in radiative J/tp + 7 decay channel 



daj/^/dy(y = 0) at Tevatron (W = 1.96 TeV) for the KMR UGDF 19], different cuts on the trans- 



verse momentum of the gluons in the loop (g§ " ) , with the on-shell gluons approximation S> qf 2 
in the hard subprocess part (denoted as "on-shell") and without it ("off-shell"). NLO correction factors 
to the hard part i^NLO UM anc ' absorptive corrections in factorized (17. 3p and in one-channel model 
(j7.ip , (j7.2j) forms are included. Contributions of all polarisations are incorporated here. 



UGDF 
approx. 


1 r ,cut\2 
\%,t 1 1 

GeV 2 


Xc(0 + ) 


Xc(l + ) 


Xc(2 + ) 


ratio 


signal 

da .]/ii-i 
£—t dy 


dcT Xc 

dy 


dvj/^-f 
dy 


dcT Xc 

dy 


dy 


d(7 Xc 

dy 


dy 


1+/0+ 


2+/0+ 


KMR, on-shell a 

(s 2 cS ) 


0.72 


41.5 


0.47 


0.51 


0.17 


0.30 


0.06 


0.36 


0.13 


0.70 


KMR, off-shell 


0.72 


10.3 


0.12 


0.06 


0.02 


0.14 


0.03 


0.17 


0.25 


0.17 


KMR, off-shell 

Seik (|7.2P, Seah. 


0.36 


29.4 


0.34 


0.13 


0.05 


0.36 


0.07 


0.15 


0.21 


0.46 


KMR, off-shell 
S cik CZ2D only 


0.36 


60.0 


0.69 


0.27 


0.10 


0.73 


0.14 


0.15 


0.21 


0.94 


HKRS [5] 


0.72 


35.0 


0.40 


0.71 


0.24 


0.45 


0.09 


0.61 


0.22 


0.73 



Q In order to avoid a confusion here, please, note this is not the rigorous on-shell case, but rather an 
approximation keeping the leading order term in q\ 2 /A^ c -expansion. 



da/dydo(y=0), nb 
1.5- 



d<j/dyd®(y=0), nb 

0.008- 




30 60 90 120 150 180 



dcj/dyd<s(y=0), nb 

0.03- 



— i , 1 , 1 0.00 

120 150 180 



FIG. 2: Distributions in relative azimuthal angle $ between outgoing protons of scalar Xc(0 + ) 
(left panel), axial-vector Xc(l + ) (middle panel) and tensor Xc(2 + ) (right panel) mesons at y = 
taking into account absorptive corrections in the factorized form (|7.3[ j (dotted lines) and in the 
one-channel model \7.1\) , {7. Sty (solid lines) forms. KMR UGDF Xm] was used. 



not ch ange d much - the effect in shapes is practically negligible. More general two-channel 
model 17|] may lead to some corrections, but they are not so big, at least, in the realistic 
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phase space regions, we are interested in 4 . 

VIII. CHARMONIUM POLARISATION EFFECTS IN THE HELICITY FRAME 

Let us start from presenting differential cross sections. As was discussed in the previous 
Section, apart from the overall normalisation of the total cross section, the absorptive cor- 
rections can noticeably change the shape of distributions in the angle $ between the final 
protons. As there is no way to measure the ^-distributions yet, the measurable effect in 
^-integrated polarisation observables, which will be discussed below, can be effectively ac- 
counted for by the gap survival factors {S^ s ) listed in Eq. ( 17. 3p . For illustration purposes, we 
show results here only for "bare" distributions (without absorptive corrections) calculated 
with KMR UGDF [19j. In principle, absorptive corrections may differ for different elements 
of the helicity density matrix, and we leave the discussion of these effects for a future study. 

In Fig. [3] we show distributions of the central exclusive x c (l + ) production cross section 
in the angle between outgoing protons $, momentum transfer t and meson rapidity y for 
different meson polarisations A = 0, ±1. In Fig. H] we present the same distributions for 
Xc(2 + ) meson production. Shapes in $ and t are rather different for both mesons - at y — 
Xc(2 + ) meson CEP is always dominated by helicity A = 0, whereas Xc(l + ) production is 
balanced by both A = and ±1 contributions, and which is the dominating one depends on 
$ and t phase space regions, we look at. 

da/dyd<J>(y=0), nb dc/dydt(y=0), nb d<j/dy, nb 




FIG. 3: Distributions of Xc(l + ) CEP (bare) cross section in relative azimuthal angle $ between 
outgoing protons (left panel), momentum transfer t along each proton line (middle panel) and 
meson rapidity y (right panel) for different meson helicities |A| = 0, 1 (dash-dotted and dashed 
lines, respectively) and for the total (summed over all X) CEP cross section (solid line). KMR 
UGDF fmj was used. 

The cross sections integrated over all possible meson rapidities y (in our case \y\ < 6.0) 
are dominated by maximal meson helicity contributions, i.e. by |A| = 1 for x c (l + ) and by 
|A| = 2 for x c (2 + ) (see Tables HT1 and II I II below) . In the last case of tensor meson the |A| = 
and |A| = 1 contributions turn out to have similar shapes in $ and t and the same order of 
magnitude, however they are quite different at y = as it is seen in Fig. HJ 

Differential distributions for fixed meson polarisations in meson rapidity y, shown in 
Figs. [3]andH] (right panels), exhibits maxima/minima in the central rapidity region y ~ 



4 We are thankful to L. Harland-Lang for helpful discussions on this point. 
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FIG. 4: Distributions of Xc(2 + ) CEP (bare) cross section in relative azimuthal angle $ between 
outgoing protons (left panel), momentum transfer t along each proton line (middle panel) and meson 
rapidity y (right panel) for different meson helicities |A| = 0, 1, 2 (short- dashed, dash-dotted and 
dashed lines, respectively) and for the total (summed over all X) CEP cross section (solid line). 
KMR UGDF 0/ was used. 



similar to one presented in Fig. HJ Interestingly enough, these maxima/minima in partial 
helicity contributions cancel each other in the total (summed over all meson helicity states) 
cross section, which has a regular and smooth behavior around y — > 0. This is due to a 
non-trivial ^/-dependence of helicity amplitudes, which is a pure kinematical effect as dis- 
cussed in Section III. Analogously, a non-trivial y-dependence of the hard g*g* — > Xe(l + , 2 + ) 
subprocess amplitudes squared shown in Fig. CE] is canceled out in sums over all helicities 
A resulting in Eq. f 13 .If) . This confirms that the appearance of non-maximal helicities is a 
kinematical effect which is absent in the spin-averaged cross-section. 



TABLE II: Non-zeroth integrated elements of the helicity matrix cr^y 1 , A = 0,±1 (in nb) for 
exclusive Xc(l + ) production at the Tevatron energy. KMR [19| and Kutak-Stasto UGDFs 18] have 
been used, absorptive corrections are not included. In the last column the ratio crjxf 1 / °~x7 l = Po^ 1 
as a measure of subleading helicity-0 contribution to the total cross section is given. 



UGDF 


rr J=1 


rr J=1 
°l.l 


rr J=1 
°1,0 


rr J=l 


Pot 1 


KMR 


0.27 


3.54 


0.51 - 0.14i 


-1.31 + 1.50z 


0.05 


KS 


0.49 


3.92 


0.85 - 0.08* 


-3.14 + 0.67i 


0.06 



Apparently, such a kinematical effect in the CEP cross section depends on cuts in meson 
rapidity in a detector. It is a result of two competing asymptotical effects: heavy (non- 
relativistic, y — > 0) meson is dominantly produced in A = state, massless (relativistic, 
y > 1) meson is produced mostly in A = ±1 states. As a result of such a competition, it 
turns out that the central exclusive x c (l + ,2 + ) production cross sections are dominated by 
the second effect such that mostly relativistic charmonia (with large rapidities y > 1) in a 
maximal helicity state are produced, so the detectors should be designed to be able to detect 
them (not only in the central rapidity region y ~ 0). 

In Tables HT1 and HTT1 we present the integrated elements of the diffractive x c (l + ,2 + ) pro- 
duction helicity matrix a{ x , = a t c P\\i- Only independent, both diagonal (purely real) and 
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TABLE III: Non-zeroth integrated elements of the helicity matrix o~ xx , 2 , A = 0,±1,±2 (in nb) for 
exclusive Xc(2 + ) production at the Tevatron energy. KMR [19|] and Kutak-Stasto UGDFs [18| have 
been used, absorptive corrections are not included. In the last column the ratio u^ 2 /a^~ 2 = Poo" 2 
as a measure of subleading helicity-0 contribution to the total cross section is given. 



UGDF 


J=2 
°00 


J=2 


J=2 
a 2,2 


J=2 
°1,Q 


J=2 
"2,0 


J=2 


J=2 
"1,-2 


J=2 
a l,2 


J=2 
a 2,-2 


Pot 2 


KMR 


0.84 


0.77 


7.58 


0.43 


0.74 


-0.54 


1.37 


1.93 


5.0 


0.05 


KS 


1.75 


1.59 


12.6 


1.02 


1.79 


-1.38 


3.17 


3.95 


8.3 


0.06 



non-diagonal (complex), non-zeroth elements are listed. Rapidity dependence of the non- 
diagonal terms (their real and imaginary parts) in the axial-vector J = 1 CEP case is 
presented for illustration in Fig. |5j 

da ,/dy, nb 



real part, X = 1 , X' = 


-1 


! imaginary part, X = 


1, V = -1 


real part, X = 1 , X' = 





imaginary part, I = 


^,x• = o 



\ !\ / 

V 

0.1 - ' v > 

o.o- *~=--vrX ...>-... 



-0.1 - 



-0.2 - 



-0.3 - 



1 1 1 1 1 1 1 1 1 1 1 1 

-6 -4 -2 2 4 6 



FIG. 5: Rapidity distributions of the independent non- diagonal elements of the helicity matrix 
daf^} /dy (real and imaginary parts) for the central exclusive axial-vector x c (l + ) meson production. 

As was already pointed out, the total CEP cross section is dominated by the maximal 
helicities. The subleading contribution of states with A = in the total CEP cross section 
is about 5-6% for both x c (l + ) and % c (2 + ) mesons. Measure of such subleading polarisation 
contribution is given by the ratio a^/a^ = Pq (see, last columns in Tables HT1 and HHP . We 

We 



calculated these ratios for two typical UGDF models: KMR [19| and Kutak-Stasto [18 
see that the ratios p 00 ~ ' are the same for both mesons, they are almost independent on 
UGDF models as the most of theoretical uncertainties (including soft rescattering effects) 
are canceled out inside them. Thus, Pq can be considered as a good model-independent 
observable for tests of underlined QCD production mechanism of \ c mesons. 

As it was demonstrated in Ref. [8j, relative contributions of different spins J are rather 
sensitive to the model of UGDF. This is due to the fact that the central exclusive production 
of Xc(l + ) is more sensitive to small gluon q t than x c (2 + ) production. Moreover, non-diagonal 
elements a^ x , , A ^ A' have much stronger sensitivity to UGDF model, than the diagonal ones. 
Therefore, experimental access to different elements of the production density matrix p J xx , (or 
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W(0,0) W(0,0) 




FIG. 6: Angular distribution W(6,<fi) of J/ip meson from radiative decay of centrally produced 
X c (l + ) (left panel) and x c (2 + ) (right panel) mesons in polar 9 and azimuthal (j) angles in the 
helicity frame. Full phase space for central exclusive Xc production is included, KMR UGDF llw 
was used. 

a \y) wou ld provide an opportunity for strong experimental constraints on the unintegrated 
gluon distributions at small x and small gluon transverse momenta qt- 

Now let us turn to the discussion of observable signal in radiative decay channel 
X c (l + ,2 + ) — > J/yj + 7. As one of the important characteristics of the production pro- 
cess pp — > piJ/iprfp, the angular distribution of the J/tp meson W(8,<j)) in the helicity 
frame is shown in Fig. [6] as a function of polar 9 and azimuthal <f> angles for both Xc(l + ; 2 + ) 
mesons. 

<W(e,<|>)> nb <W(0,(t))>^ nb 




e e 

FIG. 7: Angular distribution of J/ip meson averaged over azimuthal angle <j) as function of polar 
angle 6 (in radians) for x c (l + ) meson (solid line) and Xc(2 + ) meson (dashed line). Left panel 
corresponds to rapidity interval \y\ < 6, right panel - to \y\ < 1. 

The function W(9, 0) is a periodic one in both the angles 9 and <fi. From Eqs. ( 15 . 1 3 j) and 
(15.141) it follows that the dependence on polar angle 9 is determined mostly by the diagonal 
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terms of the production density matrix p xx , whereas ^-dependence is given by real and 
imaginary parts of non-diagonal terms. Moreover, in the angular distribution (W (9 , (p)) ^ 
averaged over cp all the terms with non-diagonal elements p xx ,, A 7^ A' are dropped out for 
both x c (l + ,2 + ) mesons. 

The ^-dependence of J/ip meson from the x c (l + ) CEP turned out to be much stronger 
than that from % c (2 + ) CEP as it is seen from Fig. O As for the ^-dependence, periods of 
oscillations for x c (l + ) and x c (2 + ) mesons are shifted by 7r/2 with respect to each other, as 
demonstrated in Fig. [7]for distribution (W(9, (p))^ averaged over (p. Also, from this figure we 
see that the amplitudes of variations in 9 become smaller and the dominance of the maximal 
helicity gets weaker and can be eliminated when one shrinks the rapidity interval in the 
phase space integral. This is a direct consequence of the specific rapidity dependence of the 
production density matrix elements in the central rapidity region y ~ 0. 



<W(9,i|))> t <W(9,(|))> < , 




FIG. 8: Angular differential distribution of J/ip meson averaged over azimuthal angle (p as function 
of polar angle 6 (in radians) and Xc meson rapidity y for x c (l + ) (left panel) and Xc(2 + ) mesons 
(right panel). KMR UGDF 0/ was used. 

Now the main question is how to measure experimentally the kinematical effects in the 
central exclusive Xc(l + ,2 + ) production reflected in the maxima/minima of rapidity distri- 
butions do~ xx /dy for different meson helicities A. Kinematical effects under consideration are 
natural in the helicity frame, in which we actually work 5 . Such a frame can be used in a 
real experiment. In principle, analogous calculations in other frames of reference can also 
be useful. As was already mentioned above, the total (summed over all meson helicities 
A) cross section is regular around y ~ and does not contain information about rapidity 
dependence of the partial helicity contributions. 

Instead, for the purpose of experimental identification of the kinematical effects and 
separation of the meson helicity contributions it is instructive to look at the correlation 
function d(W(9, (p^^/dy in both variables - polar angle of J/ip meson and rapidity y of Xc 
meson (or, equivalently, J/ip + 7 system). We present such a function for each J = 1, 2 
meson in Fig. |SJ Maxima/minima in differential distributions da xx /dy are then reflected in 
^/-dependence of the correlation function d(W(9, (p^^/dy at fixed 9. 

In principle, one can write formulae analogous to Eqs. ( 15 . 1 3[) and (I5.14p in another frame, 



5 We are thankful to S. Baranov for useful discussion of this point. 
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where the "peaks" in rapidity near y ~ can be eliminated. But then one looses some 
information about non-maximal helicity projections of the production density matrix. Since 
helicity has different meaning in different frames, our effect is purely of kinematical nature 
and thus should also be different in different frames. 



IX. CONCLUSIONS AND DISCUSSION 

Our results can be summarized as follows: 

We have calculated differential cross sections for central exclusive x c (l + ,2 + ) meson pro- 
duction for different spin polarisations in the helicity frame. The integrated cross section for 
the maximal helicity state is approximately an order of magnitude greater than that for the 
non-maximal ones. We have shown that this effect has a kinematical nature: heavy (non- 
relativistic, y — > 0) meson is dominantly produced in the A = state, massless (relativistic, 
y > 1 ) meson is produced mostly in the A = ±1 states. In the total cross section integrated 
over the meson rapidity y the second effect turns out to be dominated, i.e. relativistic with 
'j/ > 1 mesons with maximal helicity (A = ±1 for x c (l + ) and A = ±2 for x c (2 + )) are 
preferably produced in the exclusive process pp — > pXcP- 

We used two different models for the g t -dependent unintegrated gluon distributions from 
the literature - conventional KMR distribution which includes Sudakov form factor 1^ 



and the nonlinear Kutak-Stasto model [18] based on unified BFKL-DGLAP evolution. The 
contributions of the minimal helicity state A = given by the diagonal elements of the 
production density matrix Pqq 1 ' 2 are close both x c (l + ,2 + ) mesons and amount to 5 — 6 
%. They are only weekly dependent on UGDF (but rather strongly dependent on cuts 
in rapidity y). Such quantities p J Q 1 ' 2 are, therefore, good model-independent observables, 
which can be used to constrain the underlying QCD diffractive production mechanism of 
heavy quarkonia. 

We have calculated, in addition, angular distributions of J/ip meson W(6, </>) from radia- 
tive decays Xc(l + ; 2 + ) — > J/ip + 7 in the helicity frame, both analytically and numerically. 
The function W J (8,(f)) contains an important information about all independent elements 
of the production density matrix p"ly, both diagonal and non-diagonal. Azimuthal angle 
dependence of outgoing J / ip mesons is given by the real and imaginary parts of the non- 
diagonal elements p^ x ,, A 7^ A', which are very sensitive to the UGDF model used. A specific 
shape of the rapidity dependence of the J/ip meson correlation function d{W J (9,4>)) ( j ) /dy 
shown in Fig. is a result of the "peaked" shape of the Xc production density matrix 
around central rapidity region y ~ 0. Measurements of such a differential distribution 
d(W J (6,4>)) l j ) /dy averaged over azimuthal angle would allow to separate different helicity 
contributions and to identify the kinematical effects under discussion. Distributions based 
on general correlation function W J (8, <p) can, in principle, be measured at the Tevatron 
and LHC. They could also provide an independent check of the discussed diffractive QCD 
mechanism. 
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